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Abstract 

We study quantum integrable systems of interacting particles from 
the point of view, proposed in |Q]. We obtain Calogero-Moser and 
Sutherland systems as well as their relativistic-like generalization due 
to Ruijsenaars by a hamiltonian reduction of integrable systems on the 
cotangent bundles over semi-simple Lie algebras, their affine algebras 
and central extensions of loop groups respectively. The corresponding 
2d field theories form a tower of deformations. The top of this tower is 
the gauged G/G WZW model on a cylinder with inserted Wilson line 
in appropriate representation. Its degeneration yields 2d Yang-Mills 
theory, whose small radius limit is Calogero model itself. We make 
some comments about the spectra and eigenstates of the models which 
one can get from their equivalence with the field theories. 



*e-mail: gorsky@vxitep.itep.msk.su 

tpermanent address :ITEP, 117259, B. Cheremushkinskaya, Moscow, Russia 
TT~' IT' p_'pfeTirai]' nHikita@rhea.teorfys.uu.se 

1 



This paper should be considered as a natural continuation 
of 0]. The question of determination of spectra and eigenstates in 
some inegrable many-body systems will be addressed here. The main 
tool in our approach will be a hamiltonian reduction approach to the 
systems with symmetries. We shall consider two dimensional gauge 
theories, which can be treated as a topological theories and due to their 
topological nature only finite number of degrees of freedom form a 
"physical" Hilbert space of the theories. This space can be interpreted 
as a state space in some quantum mechanical problem. On the other 
hand, on can easily calculate path integrals in the topological field 
theories under consideration. It provides us with an information about 
the spectra and wave functions in the quantum mechanical problems. 
We consider first the Sutherland model [^] at some distinguished 
values of coupling constants and show that it can be obtained even at 
quantum level as a result of hamiltonian reduction from the cotangent 
bundle to the affine Lie algebra g. Then we compare this reduction 
with the reduction which is present in the two dimensional Yang-Mills 
theory and realize that these two in fact coincide. This permits us to 
get the Harish-Chandra-like formulas for the wavefunctions and for the 
kernel in the Sutherland model (actually, we use here results of Q). 
This is in the good relation with the idea of I.Cherednik to consider 
affine algebras to get exact answers for many-body integrable systems 
|11|. This result generalizes results of [M], while some methods of [Q 



could be applied in our situation too to give a satisfying treatment of 
the large A'^ behaviour of our systems. 

Recently relativistic-like Calogero-Moser-Sutherland models were 



found [12|. We show that these models (actually, their trigonometric 
version) appear in the Hamiltonian reduction procedure applied to 
the T*G (the continuation of approaches of [P,||5|), and that this re- 
duction is equivalent to that in the gauged G/G Wess-Zumino-Witten 
theory (and also in the Chern-Simons theory on a threefold which is 
a suspension of the two dimensional surface by a circle). Thus, we are 
able to tell something about the spectra and their degeneracy in these 
case. To this end we need an expression for the path integral in the 
gauged G/G WZW model on a disc. This expression can, in principle, 



be extracted from [13|. 



1 Affine Lie algebra approach to Suther- 
land model 

1.1 Hamiltonian reduction of T*g 

Let us consider some semi-simple real Lie algebra g with the Killing 
form <, > denoted also as a tr and let us fix an isomorphism between 
g and g* induced from the tr. Affine Lie algebra g is defined to be 
the central extension of the loop algebra Cg that is it is a space of the 



pairs ((/), c) where (/> is a map from the circle S^ to the g and c is a real 
number. The Lie algebra structure on g is introduced as follows: 



[{<pl{(f),Ci),{(t)2{f),C2)] = {[<Pl{f)A2{f)], / < (/>l,9i^(/>2 >), 

Now let us turn to the dual space to the g. This space g* consists of 
the pairs {A, k) where A is 5- valued one-form (actually, a distribution) 
on the circle and k is just the real number. The pairing between the 
g and g* is 

< (A, k); U, c) >= / <<t),A> +CK. 

The direct sum g (B g* we will denote as T*g. On this cotangent 
bundle the natural symplectic structure is defined: 

n= f tr{6(l)A6A)+6cA6K (1.1) 

We can define adjoint and coadjoint actions of the loop group CG 
on the g and g* . The former is defined from the commutation relations 
and the latter from the pairing, that is the element g{(p) acts as follows: 



(</.(</.), c) ^ {g{ip)cl){ip)gi^)-^ , / tr{-^g-^d^g) + c) (1.2) 

{A, k) -^ {gAg-^ + ngd^g-^ , n) (1.3) 

This action clearly preserves the symplectic structure ( |1 . 1| ) and 
thus defines a moment map 

g-* g 

which sends a 4-tuple (</>, c; A, k) to the pair {ndcj) + [A, 0], 0). 

Now let us choose an appropriate level of moment map and make 
a hamiltonian reduction under this level. To guess, what element 
of g* we should choose it is convinient to return back to the finite- 
dimensional example of this procedure. Namely, if one considers a 
reduction of the T*g with respect to the action of G by conjugation, 
then, to get a desired Calogero system, one takes an element J of g*, 
which has maximal stabilizer, different from the whole G. It is easy 
to show, that the representative of the coadjoint orbit of this element 
has the following form: 

Ju = '^ Yl (ea + e-a), (1-4) 

where u is some real number, e±a are the elements of nilpotent 
subalgebras n± C g, which correspond to the roots a, and A_|_ is the 
set of positive roots. In fact, if the corresponding Coxeter group has 



more then one orbit in the set of roots (i.e. two) then there are more 
couphng constants fll. Let us denote the coadjoint orbit of Jy by Oy 
and by the Ry the representation of G, arising upon the quantization 
of Oy. For generic J & g* let Gj will denote the stabilizer of J in the 
coadjoint orbit Oj of J, i.e. Oj = G/Gj. 

Now let us return back to the case of affine Lie algebra. We could 
just take the formula (|1.4| ) and replace everywhere usual roots by 
affine ones, but we prefer to have more geometrical reasoning. Because 
of the vanishing level k^ of the moment map, the coadjoint orbit of 
generic value of it is rather huge, while the stabilizer of the element 
(J((/5);0) E ^* is a "continious product" Hi/jeSi ^JCi/^) which is very 
small in comparison with the whole loop group CG unless for almost 
all if ^ S^ the group Gj(^p) coincides with G. 

All this argumentation shows that an appropriate value of the mo- 
ment map we are looking for is the following: 

^i={J[^llQ)■.J[^A{v) = 5{^)Jy (l5) 

Here by 5{(p) we mean a distribution ( one-form) supported at some 
marked point G S^ and normalized in such a way, that /^i 5{(p) = 1. 
The coadjoint orbit of /i is nothing, but the finite-dimensional orbit 

Oy. 

Now let us complete the reduction. To this end we should resolve 
the equation 

Kd^cj) + [A, 0] = J^{^) = 5{^)Jy (L6) 

modulo the action of stabilizer of /i, that is modulo subgroup of 
CG, consisting of g{ip) E G, such that g{Q) E Gy. We can do it as 
follows. First we use generic gauge transformation g{'^) to make A 
to be the Cartan subalgebra t d g -valued constant one-form D (it 
is always possible for non- vanishing n). We are left with the freedom 
to use the constant gauge transformations with values in the Cartan 
subgroup T C G, generated hy D - these do not touch D. Actually, 
the choice of D is not unique. The only invariant of CG action is the 
conjugacy class of monodromy exp(— D) E T C G. Let us fix some 
of these choices. Let \Xi will denote the entries of D = iX. Let us 
decompose the (^-valued function (f) on the S^ on its Cartan-valued 
part P{(p) E t and let (t>±{'^) E n± be its nilpotent- valued parts. Let 
(pa =< 4',ea >■ Then the equation (|1.6|) will take the form: 



Kd^P = 6i^)[J3]^ (L7) 

Kd^(l)a+ < D,<Pa >= Si^M]a (1-8) 

where J^ is simply Ad'^^f^JJy), [J]^ denotes the Cartan's part of J 
and 

[J]a =< J,ea >. 



From the equation ( |1.7D we deduce that D = constant and [J^]^ = 
0. This implies, that be Cartan-valued constant conjugation we can 
twist Jf, to the Jy itself. Then, ( |1.8| ) implies that locally (at if ^ Q) 
we can represent (pai^) as follows: 

0«((/j) = exp(-^ < D, a >) X M„ (1.9) 

where Mq is locally constant vector in g. Look at the right hand 
side of ( |1.8| ) leads us to the conclusion, that M^ jumps when ip goes 
through 0. This jump is equal to 

27r 
[exp( < D, a >) - 1] X M„ = [J^]„ (1.10) 

K 

The final answer for the reduction is the following: the physical 
degrees of freedom are in the exp( — f^^) and P (note that P's entries 
are purely imaginary) with the reduced symplectic structure: 

u) = —tri5P ^5X) (1.11) 

27^^ 

and we have 

exp(-^<X,a>) 



0a (V?) 



exp(-^ <X,a>)-l 
Now if we would take some simple hamiltonian system on the T*g 



whose Hamiltonian is invariant under (1^),(1^) actions, then we will 
get somewhat complicated system on the reduced symplectic manifold, 
i.e. on r*T. For example, let us take any invariant polynomial Q{(j)) 
on the Lie algebra g, that is Q G S'{g*)^ = S'{t*)'^ , where W is the 
Weyl group of G. Then we can construct a Hamiltonian on T*{g) by 
the formula 



2vr Js^ 



with somehow fixed one-form dip (for example, it will be convinient 
for us to fix a normalization ^ Jgi dip = 1). For quadratic casimir we 
get 



H2 = ^ I dip<<P,^> (1.12) 

47r Js^ 



On the reduced manifold we are left with the Hamiltonian 

//2 = -itrp2+ y ""^ (1.13) 

which coincides with the Hamiltonian of the Sutherland model. 
Here A+ denotes the set of positive roots of g. For example, for G = 



SU{N) we have the Hamiltonian for pah'-wise interaction between the 
particles with the potentiah 

yA _ 9' 



^^ sin?{xi — Xj) 
for G = SO{2N) we have: 

■' siw^yxi—Xj) sin'^[Xi + Xj) 

WA—TT^ + 



sirr{xi) sin^{xj) 

where 51,52 are couphng constants. In general case one has as 
many different coupling constants as many orbits in the root system 
has a Weyl group, i.e. two or one |^]. 

1.2 Two-dimensional Yang-Mills theory 

Let us consider some two dimensional surface S with fixed non-degenerate 
two- form LO € ri^(S) but without fixed metrics. Then one can consider 
some principal G-bundle £ over S and the space A of connections on £. 
This space is a symplectic one and group Q of gauge transformations 
acts on it preserving that symplectic structure. The moment map for 
this action in fact coincides with the curvature F = dA + A^A of the 
connection A unless our surface S has holes. In that case one has to 
specify what are boundary conditions for the gauge transformations. 

Let us take a surface S with some fixed contours Cp (in the case of 
holes one can consider non-closed contours with the ends on this holes) . 
Then for each loop C a natural gauge invariant observable is defined 
: Wc = Xa(-Pexpfp^) where a is some irreducible representation 
of G. For non-closed contours one should attach to each end of it a 
representation of G. Contours should be oriented and representations 
ar,/3r, belonging to the ends of the same contour T are to be dual. 
Then, the observable we wish to construct takes its values in the tensor 
product Qr ® af (where ar belongs to that end from this this contour 
departs due to orientation) - it is the group element P exp /p A G G in 
the representation ar, i.e. Wy = Top(Pexp/p A). This is an example 
of the coloring of the (fat) graph in terminology of Q. 

Now let us write down the action of Yang- Mills-like gauge theory 
on S. To this end we need an extra field (j) € Maps(E,g). 

Let us take any invariant polynomial function on 5, i.e. Q G 
S'{g*) . Then the action is: 

St.,q = 
= J <(t),F>+uQ{<P) (1.14) 



Let us locally decompose the tangent bundle to S to the time-like 
and space-like directions. This is a choice of real polarization - it 
has some difficulties when the surface has pants-like vertices. In the 
vicinity of generic point P E S we can always find such a coordinates 
t, X that 2-form oj will take its canonical form; 

uj = dt f\dx 

We suppose the orbits of the vector field ^ to be periodic almost 
everywhere. The symplectic structure on A (more precisely, that part 
of it, which corresponds to the vicinity of P) is 

i7 = / < 5At,5Ar > dt /\dx 



[ <6At,5A^ > 



(This symplectic structure, certainly, doesn't depend on the uj or 
whatever, it is simply fi = /j-, tr{6A A 5A).) 

Now let us suppose, that some contour T goes straight along the 
time-like direction ^. Then we get the following path integral (we 
disregard now any kinds of S except spheres with two holes 71 , 72 , 
and on each hole points Pi and P2 respectively are marked. Contour 
r goes from Pi to P2): 

j VAtVA^V^eM-Ss^y.i,Q] < vi\Peiq> j A\v2 > 

here < vi\, \v2 > are the vectors in the representations arjOp re- 
spectively. We wish to take ar to be equal R^, from the previous 
section. 

Let us remind some facts about < vi\P exp J^^ A\v2 >■ It follows 
ideologically the Kirillov constructions of the representations (and it 
was checked in [10|), that this matrix element can be represented as a 
path integral with boundary conditions. We have to take a coadjoint 
orbit Oj/ which corresponds to the R^, choose one-form 6 = d~^iJo, 
where lvq is a natural Kirillov's form on Oi, and then we can write a 
desired path integral: 

< uijPexp / A\v2 >~ / Vxexpl- {9- < At,fio >)] 

where integral in the r.h.s. needs some boundary conditions, which 
can be specified after choosing a polarization on the orbit O^, f^o '■ 
Ojj ^f g* is a, moment map corresponding to the natural G action 
on the orbit. Finally, T>x is a path integral measure, constructed 
from the 00. In the case we are interested in, (and for G = SU{N)) 
Oy = CP ~ with the symplectic form loo = Nvy. Fubini-Shtudi 
form. The corresponding moment map ^o sends the point (zq : ^i '■ 
. . . : zn-i) € CP ~^ to the matrix Jij = i'{5ij — ZiZj) G su{N)*. So, 
we can rewrite the matrtix element of the Wilson line as the following 
path integral: 



/ VziVzi exp[-z^ / ^ ZidtZi - ^ Aij{t)ziZj] (1-15) 

i ij 

Then, we put this expression into the Yang-Mills path integral and 
see, that the field At plays the role of Lagrangian multiplier, giving rise 
to the constraint ( |1.6D . Resolving it (modulo gauge transformations, 
preserving the point on the coadjoint orbit) we get the same finite- 
dimensional Hamiltonian for the eigenvalues of the monodromy around 
the space-like slice, as we had in the previous section. This proves the 
equivalence between these two theories. 

We see, that the whole construction to work, we need the orbit Oi, 
to be quantizable. This implies (in the case of G = SU{N)) that Nv 
should be integer. Further we will get more restrictive condition. 

1.3 Wavefunctions and amplitudes for the Suther- 
land model 

What can we get from this equivalence? We can solve the Sutherland 
model due to possibility of simple solving of the 2d Yang- Mills theory. 
Namely, in 2d YM theory, we can calculate exactly the path integral 
on the discomputed with the help of the form uj and a-p is the area of 
V,i.e. Jjjuj. 

In our situation we have to attach a representation Ri, to the out- 
going end of the Wilson line and the dual to the incoming one. So, the 
boundary conditions we should fix are: vectors < fi| € R*, \v2 >€z Ru 
and the monodromies 52 = -Pexp/g A, gi = Pexp/g A around the 
initial and final holes respectively. 

Gauge transformations will rotate monodromies and vectors as fol- 
lows: Qi — > hgih~^, \vi >^ Tfi^{h)\vi > i = 1,2 and the answer will 
be invariant with respect to these transformations. The measure on g^ 
is the Haar measure. Using this gauge freedom, we can make g to be 
in the Cartan subgroup T E G. The measure on T induced from the 
modding out angle variables is the product of the Haar measure on the 
T - Hi '^^i ^iid corresponding group- like VanderMonde determinant. 
To get the proper answer for the transition amplitude we should take 
a square root of this measure (to get a half- form). Finaly, note that 
in the Sutherland model itself we have a quadratic (5(</>) = |tr((/)^). 

Now let us turn to the calculation. Let us cut our cylinder along 
the contour T: we will get the disk, and on the edges, corresponding to 
r we have some group element h £ G, which should be integrated out 
with the weight < vi\Tji^^{h)\v2 >■ The whole monodromy around the 
disk is given by gihg2 h^^, so the answer for the transition amplitude 
(or evolution kernel) is the following: 

< g2; V2\ eyip[-tH]\gi; vi >= (1.16) 



From this expression we can extract the spectrum of the model as 
well as the structure of the Hilbert space (this can be achieved also 
by looking at the r*G 0): 

^ = Inv{af^ (g) a^* O R^) = 

h 

a <J)- 

h 

where <l>^ : a^ — > a^ (8) R^ is an intertwiner. 



From the answer for the kernel we can realize that more restrictive 
conditions on coupling z^ should be imposed to be possible to get the 
Sutherland system as quantum reduction: if we take gi,g2 to be T 
- valued, then the integral over h will invariant under left and right 
independent multiplication of h by the elements of T - this lead to 
the condition Tji^{T)\vi >= \vi >. It is the integral v, for which Ri, 
contains such a vector (vacuum vector of Ru)- Moreover, there exist 
only one such a vector. Let us denote it as |0 >. 

To get the expression for the wavefunction, let us look once again 
on the integral over h. It can be taken by representing character in 
the orthonormal base of representation at, so we will arrive to the 
expression: 

*a^(R}) = A/'(R})E'^™" < m\T^f^{diag[e''^])\n > (1.17) 

mn 

here J\f{{9i}) denotes the normalization factor, which comes from 
the taking into account the group-like Vandermonde determinant, tti, n 
run over the orthogonal base of the representation a and coefficients 
Cmn are defined via: 

/ dhTo,{h)mkTa{h)*nlTR^{h)QQ = CmnCkl 

so they are just Littlewood-Richardson (or Clebsh-Gordan) coeffi- 
cients. 

It is clear, that this consideration generalizes that of [^. 

2 Deformation of Calogero-Moser mod- 
els - trigonometric Ruijsenaars's system. 

2.1 Hamiltonian reduction of the cotangent 
bundle of a central extension of loop group 

In this section we will consider the deformation of the previous con- 
structions. Namely, let us consider the Hamiltonian reduction of the 



cotangent bundle to the central extended loop group G. This sym- 
plectic manifold is a space of 4-tuples 

{g:S^^G,c£ f/(l); A G n^{S'^) ®g*,K^ R). 

The action of an element h € CG is the following: 

g -^ hgh'^, A -^ hAh~^ + nhdh'^ 

K ^ K, c ^ c X S{g, h) 

where 5 is constructed from the f^(l)-valued two-cocycle on the 
group CG, T(g,h) (it provides CG with a central extension) (the co- 
cycle is supposed to be normalized : T{g,g''^) = 1): 



Sig,h) = Tih,g)Tihg,h-') 

As always, on this cotangent bundle exists a natural symplectic 
structure and the loop group action preserves it. Explicitely this form 
0, can be written in some trivialization of the cotangent bundle T*G. 
Namely, we identify it with the direct product g x G with the help 
of left-invariant one- forms kOi^ + A^^ € il^{G). The form 17 has the 
following structure: 

n= f tr[Aig-^5gf + 5A A g-'Sg] + 

I tr[nd^g ■ g-\Sg ■ g-^f - K5{d^g) ■ g-^Sg ■ g-^] + 
+c'^ScASk (2.1) 

The moment map of the action of loop group has the form: 

fi{g, c; A, k) = {gAg'^ + Kgdg'^ - A, 0) (2.2) 

It is tempting to make the Hamiltonian reduction under some ap- 
propriate level of the moment map. It is clear that in this way we de- 
form Sutherland system, presumably preserving its integrability. For 
simplicity we consider only SU{N) case. The reasoning, similar to 
that in the YM case, suggests the level of moment map to be 

nig, c; A, k) = Hj;^Id - e e+)5{(p) (2.3) 

As we did it in the previous sections, by general gauge transfor- 
mation H((p) we can make A to be constant diagonal matrix D, de- 
fined modulo affine Weyl group action, i.e. modulo permutations and 
addings of the integral diagonal matrices. Generally, such a transfor- 
mation doesn't respect the value of the moment map. Let H denotes 
the value of H(0) at the point 0. The matrix H defines a point in the 
coadjoint orbit of J, i.e. Oj = CP ~ . The semi-direct product of 
Cartan torus U{1)^~^ and Weyl group Sn acts on Oj by permuta- 
tions of homogeneous coordinates and their multiplications by phases. 
The quotient Oj/U{l)^~^ coincides with [N — 1)- simplex. 



The equation 2.2 takes the form: 



gDg~' + Kgdg-' -D = \u{^Id - / ® f+)6{ip) (2.4) 

Here / = He is some vector in C with unit norm < /, / >= 1. 
We can assume that / € R^, due to the remaining possibihty of the 
left multiphcation H — > YH with Y being diagonal unitary matrix (it 
comes from the gauge freedom which survives after diagonalization of 
A). Equation ( ^.41) yields immediately: 

g = exp(^D)G(v3) exp(-^Z)); d^G = --G5{^) 

Ki Ki r\j 

where J = 'w{j^Id — / /"*")• Let us also introduce a notation for 
the monodromy of connection D: Z = exp(— ^D) = diag{zi, . . . , zn), 
Hi Zi = 1, Zi = exp( "'^' )- We have a boundary condition: 

Q-^ZG = exp{—-)Z (2.5) 

K 

where A = e~N^ ,G = G(+0). Though the resolution of this equa- 
tion is rather interesting issue, we postpone it till the Appendix A. It 
turns out that solution has the following form: let P{z) will denote a 
characteristic polynomial of Z, i.e. 

P{Z) = \{{Z - Zi) 
i 

Let 

^±,. _ Pj^^'z) - P{z) 

^ ^ > ~ (A±^ - l)zP'{z) 

- the ratio of the finite difference and derivative of P. When A — > 1, 
the rational functions Q^{z) tend to jj. 

In these notations the matrix G can be written as follows: 

G,, = -A-V^ e''^{Q+{z,)Q-{z,)y/^ (2.6) 

A Zi — Zj 



Hgi+ij). 



sin{^) -p^ sin{ ^''^^ ) sin{- 



n 



where 9i are some undetermined phases, which are going to be the 
momenta, corresponding to the coordinates qi. It means, that on the 
reduced symplectic manifold the symplectic structure is ~ J2i dOiAdqi. 

The next step is to consider some simple Hamiltonian system on 
T*G, which is invariant with respect to the loop group adjoint action 
and make a reduction. It is obvious that any Ad-invariant function 
X : G — > R defines a Hamiltonian 



Hx= dipxio) 



For example, x±i9) = Tr(gibg^^) (here Tr is a trace in the A'"- 
dimensional fundamental representation of SU{N) ) give the Hamil- 
tonians (up to a constant factor, independent of g,) : 

H^ = Y.{e^'^±e-^'^)Y[f{q,,) (2.7) 

i j^i 

where the function f{q) is given by: 

2, . ^ . sin^{7Tiy/KN) 
^ ^^' ^ sin^nq/K) ^' 

For detailed elaboration of different limits of the Hamiltonian see, 
for example [12|. In particular, taking an appropriate scaling limit 
K — > oo will return us back to Sutherland model. Note also, that in 
order to have a real- valued Hamiltonian we need -^ < qij for any i^j, 
thus ^ should be restricted from above. 

2.2 On Gauged G/G WZW Model 



Starting from the very symplectic form (|2.lD and the moment map 
( |2.2| ) we can write a geometric bulk action for the field theory. It 
will be a sum of the term like / pdq and a moment value fixing term 
J Afndt, where Af is going to be a time-like component of a gauge 
field, serving as a Lagrange multiplier. Thus, the action is (we omit 
the term c~^dtChi, while fix the non-vanishing level k): 

S{A, g) = j dipdt tr[-A^g'^dtg - ndtg ■ g~^ ■ d^g ■ g~^ + 

+Kd-\d^g ■ g-\dg ■ g-^f) + 

+At{tig-^d^g + g"^ A^g - A^)\ 

This is the action of G/G gauged WZW model ||2^, a topological 
gauge theory recently shown to give a Verlinde formula for the dimen- 
sion of the space of conformal blocks in WZW theory on the surface 

il. 

Our theory will contain also a Wilson line in the representation Ry. 
To compute the path integral on a cylinder, as in the case of Yang- 
Mills theory,one cuts a cylinder along the Wilson line and the answer 
for the evolution cernel will have the form ( |1.16| ).In fact the represen- 
tations will run only through integrable representations of SU{N)k- 
We conjecture that the answer can be obtained also from Uq{SLi~j) 
representation theory] 18], thus establishing the relation between this 



quantum group and gauged WZW theory once again. 



2.3 Chern-Simons interpretation 

Now we present the interpretation of the previous section construction 
in terms of Chern-Simons field theory, deformed by some Hamiltonian. 
Let us consider a Chern-Simons theory with gauge group SU{N) on 
the space-time manifold being a product of the interval and a two- 
torus X = I X T^. Its action is given by: 

Scs = ^ [ tr{A AdA + -A/\AAA) (2.8) 

47r Jx 3 

We know that the phase space of Chern-Simons theory on a S x / 
with / being an interval is a moduli space of flat connections. In the 
case of inserted into the path integral Wilson lines, let us say just one 
Wilson line in the representation R, when the integral in hands looks 
like: 

VA < vi\TR{PexpJ A)\v2 > exp{-Scs{A)) 

the answer will be slightly changed, namely we will have to consider a 
moduli space of connections on the punctured surface with prescribed 
conjugacy class of monodromy of connection around the puncture. 
This class U is related to the highest weight h of the representation 
Ri, like U = exp(^^^diag(/ij)) If our surface is a two-torus S = T^, 
then we can write immediately the condition on the monodromies 
9 At 9b of connection around A and B cycles and the monordomy gc 
around the puncture : 

9A9b9a^9b^ = 9c 

Keeping in mind the condition on the conjugacy class of gc, taking 
into account that for the representation Ri, we have a sighature like u x 
diag(A^, 0, . . . , 0) we arrive to the equation ( |2.5D . The (non-rigorous) 
explanation of this relation between the connections on the two-torus 
and r*G is the following. One can consider a hamiltonian reduction of 
the symplectic space of connections on the two-torus by the subgroup 
of the group of gauge transformations, which are equal to identity 
at some (non-contractible) contour C on the torus. If one chooses 
a singular co-adjoint orbit of the moment map F = B6c, where B 
is a one-form on the C and 6c is a delta-function in the direction , 
normal to the contour C, i.e. the connection A is allowed to have a 
jump at the C, then after a reduction one is left with the limiting 
values of connection A^ and A- and with a G-valued function on 
C - monodromy g{x) {x is a coordinate on C) of connection around 
the cycle which begins and ends at the point x € C and which is 
transversal to C in homology H^(T'^,Zi). The boundary values A^ 
and A^ are related via the gauge transformation g{x), i.e. A^ = ^_|_]. 
More rigorous, but less transparent explanation involves a consid- 
eration of the space of connections on the annulus, where the gauge 



"'^We are grateful to A.Alexeiev for a discussion on this point 



group acts with a non-trivial cocycle in a Poisson brackets of its Hamil- 
tonians. This cocycle vanishes on a "diagonal" subgroup, which is a 
group of gauge transformations on the annulus, whose boundary val- 
ues are related by some diffeomorphism ip of one component of the 
boundary onto another, i.e. g{x) = g{if{x)), with x being a point on 
the one component and (p{x) is its image in the other. Then it is easy 
to see, that the hamiltonian reduction of the space of connections by 
the action of this group can be done in two steps; first, we reduce by 
gauge transformations, equal to identity at the boundary of the annu- 
lus, then we take a symplectic quotient with respect to this "diagonal" 
subgroup, isomorphic to the loop group CG^. 

Also we could point out the equivalence of the sector of the space 
of observables in the Chern-Simons theory on S x S*^ , namely, observ- 
ables, which depend on At with dt being a direction, tangent to S^ 
and the gauged G/G WZW theory. This can be proved either consid- 
ering a gauge fixing At = const^ and diagonal, or by poining out that 
the wavefunctions in the Chern-Simons theory are those in the theory 
with an action J dt /^ tr^AdtA) projected onto the space of function- 
al, invariant under the gauge group action. This involves taking into 
account the non-invariance of the polarization jj and this is why the 
WZW part of the action appears (see [|l3| for detailed discussion). 

2.4 The spectrum and wavefunctions of the 
model 

From this equivalence we can derive a spectrum of the Ruijsenaars 
model. Namely we will show that it is contained in the spectrum of 
the free relativistic fermionic system on the circle. By free relativistic 
model of fermions on the circle we mean the system of particles on the 
circle with coordinates < qi < 1 and conjugate momenta (rapidities) 
6i living also on the circle, with the Hamiltonian 

i 

and symplectic structure 

zvr 

(in fact, the shift k ^ k -|- A^ is a quantum effect and can be deduced 
from [O], so we don't consider its derivation here). The fermionic na- 
ture of these particles amounts to the phase space (T^"-*^ x T^'^)/Sn 
- just the moduli space of flat connections on the torus. Here T^~^ 
is a Cartan torus for SU{N) (center of mass is fixed) and Sn is a 
symmetric group. We call them fermions, because their wavefunction 
vanishes on the diagonals Zi = Zj in the N — 1-torus (it contains a 
factor like A(e'^*) - a group-like Vandermonde) . It is obvious, that 



We are grateful to V.Fock for valuable discussions on this point 



the corresponding Schrodinger equation will be a difference equation, 
because cos{-^^dq) is a finite difference operator. Its eigenfunctions 
on the circle are exponents e^'^'"''^ (n is necessarily integral) with eigen- 
value 

27rn , 

and the total spectrum will be given by the sum 

i 

where we have to impose the following conditions, coming from the 
invariance En = En+K+N and the symmetric group action: 

(k + A^) > nAf > . . . > rij > . . . > ni > (2.9) 

Y,ni<{K + N) (2.10) 

i 

In our situation fermions on the circle are not completely " free" and 
the spectrum is restricted. For our purposes will be more convinient 
to form a vector /i = (nj — z), i = 1, 2, . . . , A^. It will correspond to the 
dominant weights of the SU{N) irreducible representations. Now we 
can write an answer for the spectrum: 

i?.=Vco.(^^) 
^ ^ K + iV^ 

and h satisfy the selection rule: there exists non-trivial intertwiner 

^h : «/» ^ «/» ® ^^ 

(essentually it implies, that h = A + i^p, where Aj > 0, i = 1, . . . , A^. 

An explanation of this spectrum can be derived from the properties 
of gauged GjG WT^N model |l^. Namely, to get a partition function 
of our interest we have to compute a path integral on the two-torus: 

fvAVgexp{SGWZw{9,A) + / Mg + g-^))xRAP<^^P [ ^) 

(2.11) 
where C is a non-contractible contour on the torus, XRi, is a char- 
acter in the representation R,^. By the standard arguments we can 
calculate it by cutting the torus along C and computing first the path 
integral in the theory Scwzwig, A) + /Tr(5f + g~^) on the annulus 
with coinciding monodromies h of the restrictions of the connection 
A on the boundaries of the annulus and then integrate the answer 
over H with the weight xr^, (h)- To make a further progress one needs 
some information about the path integral on the annulus (or cylinder, 
what is the same). In principle, this can be deduced from |jl3[. We 
will not present here the whole derivation, just note that the path 
integral on the cylinder can be localized (either by some equivariant 



cohomology-like technique for SU{N) gauge group itself or by using 
a gauge, where g is diagonal and then use the usual localization for 
abelian gauge theory |]l^]) on the g with special values of eigenvalues, 
namely those which correspond to the integrable representations of 
SU{N)k, i.e. with the signature, given by numbers n^, satisfying the 
same restrictions (p.g|),( p^ . In arbitrary gauge theory in two dimen- 
sions the path integral on the disk and on the cylinder can be written 
in terms of sum over all irreducible representations of the gauge group 
via the characters in these representations. We conjecture that tak- 
ing into account Wilson line in the representation Ri, will lead us 



to Macdonald polynomials |18]. We hope to return to more detailed 
investigation in the future. 

3 Conclusions and discussion 

In this paper we have considered some particular integrable defor- 
mations of Calogero-Moser systems, which preserve their trigonomet- 
ric form, but make the quantum Schrodinger operator a finite differ- 
ence operator. We have observed an equivalence of this theory with a 
gauged SU{N)/SU{N) WZW model with Hamiltonian on the cylinder 
with inserted Wilson line. The level of WZW model and the repre- 
sentation in which Wilson line is taken are encoded in the coupling 
constant of the model and in the mass scale. 

There are some important and intriguing questions, which remained 
beyond the scope of this paper. The first one is the elliptic analogue of 
the topological theories like gauged WZW or its degeneration - Yang- 
Mills theory. We expect this elliptic generalization to exist, because 
in the corresponding quantum integrable systems of particles there is 
such an analogud^. When elliptic curve degenerates one gets trigono- 
metric version of integrable model. One of the possible ways is to 
consider double-loop algebras [19|,[17| and/or Sklyanin algebra and 



its multivariable generalizations as the algebras of symmetries of cor- 
responding theories. Some indications of the appearence of Sklyanin 
algebra are in the Askey- Wilson polynomials, which are particular ex- 
amples of the wave- functions of systems under consideration p^ ,|14|. 



Another interesting question concerns investigation of different 
limits of our theories when coupling constants go to some extremal 
values. In fact, we have two important parameters in hands: 



and 



2iTi 

QK + N 



t = g'^+i 



We have conjectured that in the gauged G/G WZW theory with 
the Wilson line in the representation R^, the answers for path integrals 



^See pO| for the progress in this direction 



on the surfaces with boundaries can be expressed in terms of Macdon- 
ald polynomials Pf^{x;q,t). These polynomials have very interesting 
limits (for example, they are responsible for the zonal spherical func- 
tions on the p -adic groups, namely Mautner-Cartier polynomials when 
q = 0,t = 1/p) if q, t tend to some exceptional values. Unfortunately, 
in our approach q should be a root of unity and t should be an integral 
power of q, so there is no obvious way how to recognize all this beauty 
of Macdonald polynomials in gauged G/G WZW theory for SU{N) 
group. Presumably SL{N, C) can overcome this problem. 

It is known that Ruijsenaars system is closely related with the 
solitonic sector of Toda like theories, in particular asymptotics of the 
wave funstions provides S-matrix for scattering in the solitonic sector 
|23|. The qualitative explanation of this relation looks as follows. La- 
grangian approach for Toda theories is based on GWZW actions with 
some additional dependence on the spectral parameter. In the consid- 
eration above we added Hamiltonians and choosed some peculiar level 
of the momentum map. Thus the final action is nothing but the gen- 
erating functional for GWZW while the level of the momentum map 
fixes the particular solitonic solutions. Hence we considered in essence 
the generating functional and solitonic S-matrix naturally appeares in 
this context. Certainly this relation needs for further clarification, es- 
pecially in context of (q)KZ equations describing the soliton scattering 
amplitudes. We plan to consider this problem elsewhere. 

Recently, in the works of Dunkl, Opdam and Heckmann ||] a pow- 
erful method of solving Sutherland and Calogero - type models was 
found. It is based on the introducing of operators acting on the group 
algebra of the symmetric group C[Sn]- This method was applied also 
in a seria of papers to a more general systems, including spin long- 
range interactions (see, e.g. 0). To obtain a meaning of these opera- 
tors in a gauge theory remains to be quite important question. They 
could be related somehow to the problem of Gribov copies in a gauge 
fixing in a gauge theory and, therefore can be of some improtance in 
a four-dimensional gauge theory. 

Another important direction is the finite-dimensional verification 
of equivalences, obtained between finite-dimensional system and (at 
first sight) infinite-dimensional. In the case of Sutherland model one 
can check, that there is an intermediate step in the procedure of Hamil- 
tonian reduction, leaving us with the cotangent bundle to the group 
SU{N). In the case of afhne group instead of affine Lie algebra the 
situation seems to be slightly complicated and one hopes to get as 
a finite-dimensional analogue of T*SU{N) one of the SU{N) dou- 
bles, most probably the Heisenberg double, which has an open dense 
symplectic leaf. But, instead of Hamiltonian action of SU{N) on 
its cotangent bundle in this situation we have a Poissonian action of 
Poisson-Lie group, which is slightly more involved situation Q, |jl6|. 
Still it seems to be interesting question to answer. We think that 
in this way one can understand the relation of the quantum group 



representation theory to our problems |18]. 



Let us finally note that finite-dimensional Ruijsenaars system can 
be consistently quantized. To this end one should find a proper nor- 
mal ordered form od the Hamiltonian.For the N = 2 case it has the 
following form 

H = $+T_$_ + $_T+$+ 

where $+ = $1, $_ = $2-In our approuch we have the very normal 
ordering which preserves the integrability. 
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Appendix A 



In this section we prove the relation (|2.q ). 

The necessary condition for the equation ( |2.5| ) to hold is that the 
eigenvalues of matrices Z and 



exp( J)Z = \{Id + (A"'^ - 1)/ ® r)Z 



coincide. 

It leads to the system of equations 

N 



^i 



k = l,...,N; <!>, = I/i|2 

Also we have to satisfy the normalization condition J2i ^i = 1- Let 
us note that our equations ( |0.1| ) resemble those arising in the defini- 



tion of the Gelfand-Tzetlin coordinates |15]. The reason is obvious: 
Gelfand-Tzetlin coordinates are related to the imbeddings of S{U{N — 
1) X [/(!)) into SU{N), the quotient SU{N)/S{U{N - 1) x ^7(1)) is 
precisely CP ~ - our coadjoint orbit. Also let us note that the so- 
lution of these equation (written on a next line) looks like expression 
in the right hand side of the Bethe ansatz equations - the reason for 
this is to be understood. Our claim is: <I>j = Q~^{zi). We prove this 
formula by twofold evaluation of the following contour integrals in the 
complex plane around infinity: 

1 / ,,.f(A.) 



and 



2m Jcac zP{z) 



' ^ d=- ^<^^' 



2»i 7c„ {z - A-lz4)-P(--) 
Contour Coo is oriented counterclockwise and goes around infinity. 
Both of them are obviously equal to A (residue at infinity). On the 
other hand, Iq^Iu are given by the sum of residues at poles of the 
integrands inside of Coo • It gives 

/o = l + 5:(A^-l)Q+(z,) 

i 

A^-1 



4 = EQ^(^i)f 



z^x-i 



and we obtain j This is equivalent to X{Zvi + (A —!)</, Zvi > 
f) = ZiVi, which leads to 



N iN/^ x-1 v\-l 



Vi = gi{X-'' -l){z,X-'-Z)-'f 



where gi are some constants to be found from the normahzation 

<Vi,Vj >= 5ij 
To this end we have to calculate an expression 



^r' = \9i\-^ 



(2_A^-A-^)^ 



^-i 



{z,\-^-z,){zr'\-zT') 



—^ (0.2) 



This expression can be evaluated with help of residues of the fol- 
lowing contour integral 



J.^^ 



P{Xz) 



27ri Jc^ {z - z,\-^YP{z) 

It is obvious that J^ = 0, because at infinity the integrand behaves 
as z^^ . On the other hand, it implies: 



P{\-^Zk) 



Y^ Res,=^^ 



P{\z) 



[z-zi\~^YP{zy 



(0.3) 



thus, 



Now we are in a position to write down an expression for matrix 



iv-i A' 



'N 



Zi — X~^z 



^3 



G 



Finally, after substituting here (^,(^) we just get (|2J). The 
phases 9i are restricted according to the condition det{G) = 1, so this 
is a condition, imposed on the X^j ^i; which doesn't take us here. Still, 
for completeness, let us calculate this determinant: 

det{G)= E(-rnG'.,<.« = 



e^S/^AVn(Q+(z,)Q-(z,))^/^ E (-rllT^ 



A-^-1 



<76Sjv i ^ ^* ^'^(*) 



e*E,«.A^(A-^-l)^ 



(A-l)(A-^-l) 
(A^-l)(A-^-l) 



nl/2 



\[{Q+{z^)Q-{z^)f/^det\\—^^ 



Zi — A Zj 



(Cauchy identity) 



e '^3 •'A 2 



-N 



1 



'^n-nzi^r.n 



{Zj - Zkf 



) ;<fc(^J' ~^k^ ^)'y^3^ ^ -^k) 



X 



{Xzj - Zk)iX ^Zj - Zk) 



(Zj - ZkY 

= e ^j ^ = 1 

The determinant det\\ ,_}-i , || can be evaluated by representing it as 
a correlator of free fermions on a complex plane inserted into points Zi 
and Xzi - it is a "physical" proof of Cauchy identities used also in |T^^. 



So, the answer for this determinant cancels precisely all Zi dependent 
factors in the product (this could be also deduced from the absence of 
poles of this meromorphic function, therefore, it should be constant). 



Appendix B 



Here we derive an expression for the wavefunctions of Calogero 
model via the large k limit in those for Sutherland model. We start 
from the expression for the transition kernel 

T 
}C{g2,V2\gi,vi) = ^exp( 2 ^ ("))^" (5^2, 1^2151, ^'l) 

a ^ 



K^a{g2,V2\gi,vi) = dhxoiihgih ^g^^) <V2\TR^{h)\vi> (0.1) 

JSU(N) 
Let us write 51^2 = exp(^Qi 2), then hgih^^g2^ = exp(^<l>) with 

^ = hQih-^-Q2 + 0{-) 

K 

Now we can write down an expression for the eigenvalues of the 

matrix hgih^^g2 as e~^ where ^j are eigenvalues of <1>. By the Weyl 
character formula we have: 



dim{a) 



A{i) 
detWPfl 



X.iU) = -^^ (0.2) 

where U is conjugate to diag{(3i, . . . , I3n)- Here collection of pi's is the 
highest weight h of the irreps a^. 

The formula we will derive uses the representation of the large 
K limit of the normalized character through the Harish-Chandra(- 
Itsykson-Zuber) integral, namely, if - = P is fixed, then 

lim ^^^= I dheMiTr(PhQh-^) 

K^oo dim{a) Jsu{N) 



We thank A.Marshakov for pointing out this useful application of Wick theorem 



where exp(z^) = diag{U), and diag{U) denotes any diagonalized form 
of the matrix U . Substituting this expression into the integral in ( |0.lD , 
we get an integral of the type: 



/ / dhdh' exp{iTr{h'{hQih-^ - Q2)h''^P) 

J JSU{N) 

<V2\TR,Xh)\vi> (0.3) 

Finally, we can extract a wavefunction from the condition: 

^%{Q2)-^p{Qi) =11 dhdh'}Cp{glv^\g^',v1') 

^ J JSU{N) 

where g'^ = hgh~^, v'^ = Tji^{h)v. 
It yields essentually : 

A(P)A(Q) / dheicp{itr{PhQh-^) < 0\TR^{h)\0 > 

JSU(N) 

In the case of A^ = 2 this integral reduces to the integral represen- 
tation for the Bessel function. It is interesting to note, that the initial 
expression of the wavefunction for the Sutherland model could be re- 
garded as the same formula applied to the central extended loop group 
- in that case we would get some path integral, which we can evaluate 
using, for example, localization technique, or, by choosing appropriate 
coordinates on the co-adjoint orbit of affine algebra ( [p^] ). 



